Numerical analysis of the existence and stability of nonlinear excitations in a 
parametric model of ferromagnetic chain. 
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A parametrized spin model was recently introduced and intended for one- dimensional ferromag- 
nets with a deformable Zeeman energy. This model is revisited and given more realistic interpre- 
tation in terms of a model for ferromagnetic systems with nonconvex anisotropics. A main virtue 
of the improved form is its exact reduction to the discrete Remoissenet-Peyrard model, i.e. a 
parametrized version of the Takeno-Homma's discrete sine-Gordon model. The spin-wave phase of 
the improved parametrized spin model is investigated assuming both harmonic and anharmonic exci- 
tations. Intrinsic-self-localized modes, regarded as zone-boundary breather spin waves, are pointed 
out by simulating the nonlinear difference equations describing the spin equilibrium positions in 
the chain, and are shown to exist irrespective of values of the model parameter. Domain-wall tex- 
tures of the model are also numerically examined in terms of kink solitons and with regard to the 
parametrization. 

PACS numbers: 61.72.Bb, 62.30.-fd, 63.20.Ry. 



I. INTRODUCTION 

One-dimensional(lI?) magnets are simplest in their 
description while displaying at relatively low magnetic 
fields, a continuous symmetry related to spin preces- 
sions with respect to the chain axis [J. These sys- 
tems have been widely investigated over the past years 
to both viewpoints of large and small-amplitude spin 
displacements and in connection with several real mag- 
netic compounds, two most considered been TMMC and 
CsNiFs 2, 3., _4, ^_5.,, 6J. However, one-dimensional mag- 
nets are rather aboundant in the nature therefore imply- 
ing quite various and often complex magnetic features 
e.g. structures. The first, CsNiF^, is probably the best 
studied and relatively clear understanding of its struc- 
tural properties has now been gained. In this compound, 
the Nickel ions align into a linear chain and are cou- 
pled through their magnetic moments giving rise to a 
ID ferromagnetic system. In connection with this par- 
ticular arrangement, magnetic moments experience an 
anisotropy field so called single-ion anisotropy and due 
to the crystalline structure. The crystalline anisotropy 
favors in-plane spin precessions where the rotation plane 
can be perpendicular to the chain axis. 
When no magnetic field acts on the system, spins can 
freely rotate and their precessions lead either to large- 
amplitude or to small-amplitude exictations. Large pre- 
cessions imply complete rotation hence 27r excitations 
which in theory are described as solitons. Solitons re- 
late to domain walls or the elementary excitations of the 
collective spins. To a physical viewpoint, single domain 
walls represent local magnetic orders which, under ther- 
mal fiuctuations, can nucleate giving rise to either short 
or long-range ordered ferromagnetic phase depending on 
the strength of the thermal fluctuations. As opposed to 
domain walls, small-amplitude magnetic excitations are 
more extended and dispersed hence do not necessarily 
cause complete breaking of a prevailing magnetic order. 



At the most they will promote order-disorder instabili- 
ties governed by spin-wave dispersions. However, at low 
enough temperatures, spin waves can coherently order 
into anharmonic or breather-like excitations. In fact, 
these breather-like structures form an hierarchy of varied 
localized modes in the ferromagnetic groundstate 7] to 
which belong the intrinsic-self-localized(ISL) modes re- 
cently discussed H, 0, 0| for both ferromagnetic and an- 
tiferromagnetic compounds. Since they describe excited 
states of the ferromagnetic groundstate, anharmonic lo- 
calized modes can be looked on as classical anharmonic 
fiuctuations in this groundstate. In general, such anhar- 
monic fluctuations will exist irrespective of the wavector 
but depending mainly on the magnitude of the spin dis- 
placement. However, ISL modes are particular in that 
they need the wavector be exactly at the edge of the 
first Brillouin zone such that neighbour spins coherently 
oscillate anti-phase them and others but at the same fre- 
quency. In fact, all these coherent anharmonic objects 
extend our eyesight of the spin-wave phase of II? mag- 
nets beyond the usual one restricting to linear phenom- 
ena. 

In the presence of a uniform field perpendicular to the 
chain axis, the spin rotations become strictly planar and 
the ferromagnetic order almost confined along the chain 
axis. As a consequence the spin equilibrium and dynam- 
ical properties will strongly dependent on the magnitude 
of the applied field. While solitons in this last context 
still correspond to the progressive 27r rotation of spins 
about the easy axis, the region spanned by the spins is 
now determined by the field. Refering to the domain-wall 
texture, this region corresponds exactly to the spatial ex- 
tension of a single domain hence the soliton characteristic 
width. Analytically, one finds that this region varies as 
the inverse square root of the magnitude of the applied 
field. Moreover, for the discrete planar model largely 
discussed in the literature, the magnetic field determines 
conditions for the propagation of soliton along the easy 
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axis via a defined dispersion law as Takeno et al. pointed 
out [Hill- 
Very recently, a ID parametrized spin model appeared 
in the literature Ts!, and was shown to reduce to 
a parametrized sine-Gordon mo del p reviously suggested 
by Remoissenet and Pevrard [iM. [l7l | and studied by sev- 
eral authors [l^ IT^ I2(tI I21I l23|T In connection with this 
well-known model, the nonlinear dynamics of the ID 
parametrized spin chain has been investigated 0, 
and reveal rich topological soliton features. In spite of 
a lacke of sound experimental contexts conforting such 
kind of parametrization, the existence of exact soliton 
solutions is a relevant source of motivation for an inter- 
est to this new model. Keep in mind that none of the 
current models i.e., neither the isotropic nor the planar 
ising or Heisenberg models has yet demonstrated ability 
to give complete account of experiments on the available 
magnetic systems, and will certainly not be more help- 
ful for an understanding of the physics behind structures 
shown by the various magnetic alloys. In its original 
form, the parametrized spin model describes a ID fer- 
romagnetic chain subjected to an easy-plane single-ion 
anisitropy and a Zeeman energy with nonconvex depen- 
dence on the spin variable. While the author [3| relates 
this nonconvex dependence to some deformability, a more 
efficient and effective interpretation can be given in terms 
of a non-conventional crystalline structure such that the 
magnetic field cannot couple exactly to none of the spin 
components. In this context, the parametrization can 
stand as a relevant theoretical mean to include devia- 
tions from conventional crystalline structures. We there- 
fore suspect such an approach to provide an interesting 
way to model some real complex structural magnets ob- 
tained e.g. by isotopic substitutions, magnetic dopings, 
alloyings, etc. While agreing with the original modelling, 
we retract from the explicit inclusion of the cristalline 
anisotropy energy. Indeed, the nonconvex Zeeman term 
is to our view already self-sufficient i.e. involves energies 
from the two main anisotropics. 

We will focus on one single point i.e. the static prop- 
erties of spins in this model, the ultimate goal being to 
check the feasability of such a parametrized spin model. 
In particular we will probe possible analytical weakness, 
namely whether some restrictions exist that prevent har- 
monic as well as anharmonic spin-wave and kink soliton 
phases. 



II. THE SPIN- WAVE PHASE. 

Consider a 1_D magnetic system with ferromagnetic in- 
teractions among the spins. We assume that the mag- 
netic field does not couples exactly to none of the spin 
components such that for a field directed perpendicular 
to the z axis, the system Hamiltonian writes: 



where J is the exchange coupling and r some real param- 
eter. h(r) is a function defined as: 

^2 



h{r)^hoil~ry 



(2) 



where ho is the effective uniform field that, for a well- 
defined anisotropic crystalline structure with easy axis z, 
should act on the component of the spin vector s„. 
Compared with the original model > the new form 
retains only the nearest-neighbour exchange energy and 
the nonconvex anisotropy energy. The reason for such 
simplification follows from the remark that one can ex- 
pand the nonconvex energy for weak values of r(or more 
precisely, of the ratio 2r/(l -1-7"^)) and up to the term in 
(s^)^, both Zeeman-like and single- ion-like contributions 
are generated. So in its present form, the parametrized 
model is self-consistent at least for the two main encoun- 
tered anisotropics in ID ferromagnets. Moreover, fur- 
ther expansions give new energy terms each either con- 
tributing of a new symmetry or enforcing an already ex- 
isting symmetry. To this viewpoint, such an expansion 
gives rise to a generalized free-energy-like Hamiltonian 
in which the first term(i.e. the exchange energy) can be 
looked on as the distorsion energy. 

Because of the product of the parameter ho with all terms 
of this expansion, a relevant question arises as to how 
symmetries get broken in the system. The answer re- 
sides in the fact that coefficients of this expansion are 
not of the same power in r. Therefore one can always de- 
fine local fields characterized by different powers of the 
deformability parameter r, which then appears to rule 
symmetry-breaking phenomena in the non-conventional 
crystalline structure. In view of all these remarks, we 
will avoid calling the last term in |^ a Zeeman energy 
as in the previous works. Instead, we will refer to as a 
nonconvex anisotropy energy. 

As we are interested with the spin-wave phase of this 
model, we must distinguish two classes of spin waves 
characterized by two distinct amplitude regimes. How- 
ever, both spin-wave regimes lie in the ferromagnetic 
groundstate, of which they describe classical fiuctua- 
tions. One of these two regimes is dominated by small- 
amplitude deviations from the uniformly ordered states 
Sn = S, i.e. involves harmonic waves. This regime is 
best described using coherent state variables as given by 
the Holstein-Primakov transformation j^^. Expressing 
the spin variables s„ as functions of these coherent-state 
variables, next expanding the nonconvex term in r and 
grouping all quadratic terms in the coherent states we de- 
rive the following dispersion spectrum for the harmonic 
excitations: 



uj{q) ~ ujosin (k a/2) + — 
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H = -2J^ s„s„+i + h{r) - 
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Where a is the lattice spacing and k the wavector. This 
dispersion relation reveals a gap in the groundstate of the 
spin- wave phase i.e. at fc = 0, of amplitude: 

2 

(4) 



1 



3 



Let us expand the function V{s^,r) as a polynomial in 



FIG. 1: The spin-wave dispersion spectrum, plotted in the 
reduced wavector ka/ir. 



and which closes continuously as r increases. Accord- 
ing to figure^ the field-induced gap causes uniform shift 
of the whole excitation spectrum such that not only the 
fc = spin states of the ferromagnetic groundstate will 
be affected by its closing. Still, excitations involving 
the spin states near fc = will be more prone to har- 
monic phenomena as opposed to excitations of the spin 
states lying close to the edge of the first BriuUoiun zone, 
from which relatively large-amplitude phenomena are ex- 
pected. Since these last excitations suggest possible an- 
harmonic spin waves the classical spin motion should 
obey discrete nonlinear equations. Owing to the com- 
plex form of the model Hamiltonian, the technique of 
derivation of these equations deserve attention. We will 
follow a previous approach transforming the Hamilto- 
nian as a function of scale (rising and lowering) variables 
= ± «s^, so that can be rewritten: 



H 
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1 + + 2rsl 



(5) 
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The next step is to apply the following canonical equa- 
tion: 



which leads to the formal equation: 



(7) 
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dt 



= -2J {s+^, + - 4 «+i + 

+ [stVis^^r)] . (8) 
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2r 



1 + r2 



(9) 



Proceeding term after term, we can establish that the 
commutator of the rising variable with the polynomial 
function Vi of ultimately reduces to: 



[4,yi«,r)] = [4,< 



dVi 
dsl. ' 



(10) 



With help of these transformations, the discrete equa- 
tion (jHl becomes more explicit. However, to derive 
this explicit form s+ must first be expressed as a time- 
dependent quantity. In this goal we define: 



s„(t)e 



i(kna — Q,t) 



(11) 



Inserting Hll|l in JSJ and (|10|l . and next separating real 
and imaginary parts of the resulting equation, we arrive 
at two new equations i.e.: 



dT 



= -Vl - |s,i(t)|2 [sn+i{t) - Sn-i(t)] sin(fc). 



2J 



(12) 



v/l-|s„+i|2 + v/l-|s„-l|2 



\/l - |s„|2 [Sn+1 + Sn-l] COs{k) 

h{r) s„ 



2J 

h 



1 + r2 -f 2rVl- |s„|2 



2J ' 



(13) 



where we assume S = 1. Equations (|12|l and (|13|) form 
the complete set governing the dynamics of anharmonic 
spin-wave excitations of the parametrized model. How- 
ever, to keep our primary goal we will concentrate on the 
static part of this set. Recall that this static part is most 
relevant for the equilibrium properties of the ferromag- 
netic groundstate in the presence of the anharmonic clas- 
sical(spatial) fiuctuations. In this context, the discrete 
static equations H13|) provide reccurence relations be- 
tween positions of neighbour spins fluctuating about their 
equilibra. The presence of O in these equations, which is 
the frequency at which spins fluctuate oscillating in time, 
traduces an account of adiabatic spin fluctuations thus 
meaning quite rich equilibrium properties. Notably, ISL 
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modes will occur as anharmonic classical fluctuations cor- 
responding to the spin- wave phase in which Q gets closer 
and closer to the cut-off frequency uj{k — vr/a), such 
that spins of same partity(e.g. even-site spins) coher- 
ently oscillate in phase being completely decoupled from 
the spins of opposite parity(odd-site spins). So to say we 
expect two co-existing anharmonic modes at this specific 
wavector, leading to two degenerate localized spin-wave 
excitations in the ferromagnetic groundstate. Here the 
degeneracy relates to the fact that both localized spin- 
wave modes fluctuate at a common frequency though be- 
ing completely decoupled one from the other. Instruc- 
tively, the phenomenon just described looks quite like the 
nonlinear self-localized objects predicted by Chen and 
Mills '23| in photonic band-gap crystals, and extended 
to other lattice systems e.g. atomic systems (25t . i2& i27J | 
and nonlinear transmission lines |23, 123 . The most at- 
tracting feature in these particular anharmonic modes is 
their common breather-like shape profiles. Stress that 
by breather we understand an harmonic wave hidden in 
a spatial nonlinear topological shell. Being robust by 
virtue of its topological feature, we readily expect this 
spatial modulation to increase lifetime of the harmonic 
wave, namely by fixing the frequency of its "breath" in- 
side the shell. In this respect, the most common be- 
lief supported by numerous theoretical results is that the 
breather frequency will depend on the width of the spa- 
tial topological shell, which often is one of the two most 
popular topological solitons i.e. kink or pulse. 
Let us now put the previous discussions on a firm ground 
by carrying out numerical simulations of the discrete non- 
linear equation (|13|l . As a starting note, following these 
previous discussions we can also term ISL modes zone- 
boundary breathers. Indeed such a description reflects 
rather well both their well-known waveshapes, and the 
fact that they involve wavectors lying close to the bound- 
ary of the first Brillouin zone(fc = n/a). In connection 
with this last requirement, it is useful to clearly define 
the interval of allowed frequencies as given by the disper- 
sion relation From the condition < k < n/a the 
dispersion relation yields: 

h[^y<n<.^ + k{l^)\ (14) 

As it turns out, placing ourselves near the edge of the 
first Brillouin zone means choosing Q close to the up- 
per limit of the constraint H14() . Thus, once we fix the 
ratio ho/ J an appropriate choice of the zone-boundary 
breather frequency follows. In our simulations, ho/ J is 
kept fixed at 0.1 which allows us selecting a frequency 
of about 4(J = fi = 1) or slightly less. Last, periodic 
boundary conditions will be assumed and applied such 
that neighbouring spins start fluctuating simultaneously 
with opposite amplitudes. Numerically, this amounts to 
set Sq'' = — Sq^ where superscripts "ep" and "op" refer to 
even-parity and odd-parity modes, respectively. 
On flgures |21 and 13 we illustrate shapes of the odd and 
even parity modes for the conventional planar model 
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FIG. 2: Odd-parity zone-boundary breather: ideal planar 
model r = 0. 
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FIG. 3: Even-parity zone-boundary breather: idela planar 
model r = 0. 

r = 0. We clearly see that their maxima are uniformly 
shifted which is a signature of their anti-phase oscilla- 
tions. Otherwise, the twin character of their shapes con- 
stitutes a sound proof of a ferromagnetic ordering in the 
groundstate of the conventional planar model. An in- 
structive step following this previous relevant one is to 
assume weak values of r and check the becoming of the 
zone-boundary breathers in flgure [21 and |3| Figures 01 
and El show numerical results for r — 0.05. As these 
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FIG. 4: Odd-parity zone-boundary breather: r = 0.05. 
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FIG. 5: Even-parity zone-boundary breather: r = 0.05. 
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FIG. 6: Shape of the zone-boundary breather: r = 0.2. 



figures suggest, a nonzero value of r does not clianges 
the physics contained in the conventional planar model. 
Namely, the groundstate remains ferromagnetically or- 
dered and manifestly as stable as the groundstate of 
the conventional planar model. In fact, r only affects 
shape profiles of the odd and even parity zone-boundary 
breathers but not their existence. More explicitely, an 
increase of r causes their spatial extensions (widths) to 
increase while preserving their twin structures. From 
this step we can already conclude that the parametrized 
model allow smooth variation of characteristic magnetic 
features from the ideal limit corresponding to the usual 
planar model. However, to permit global quantitative 
and qualitative appreciations we will display numerical 
solutions for larger values of r. Thus, figures El [3 
and llUl are some representative shape profiles. As we are 
aware of the twin character of the two modes, only odd- 
parity solutions are shown in this last set of figures. 



III. TEXTURES OF DOMAIN WALLS. 

Soliton solutions of the Remoissenet-Peyrard [T^ [TtI I 
model, which we recall is currently the most celebrated 
generalized sine-Gordon model, have been thoroughly 
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FIG. 7: Shape of the zone-boundary breather: r = 0.5. 
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FIG. 8: Shape of the zone-boundary breather: r = 0.8. 

investigated [11 El S El 113 both analytically and 
numerically. Various applications have so far been en- 
visaged extending from Biophysics to crystal growths 
through magnetic excitations in II? magnets. The dis- 
crete version of the sine-Gordon model has been intro- 
duced two decades ago 11|, iljj as a model for DNA con- 
formational transitions. This model was next extended 
to various other contexts but the most evident system in 
which it naturally arises is the ID anisotropic Heisenberg 
model. Here we will see that the parametrized model ^ 
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FIG. 9: Shape of the zone-boundary breather: r — 1. 
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FIG. 10; Shape of the zone-boundary breather: r = 10. 



FIG. 11; Kink lattice of the planar Heisenberg model(r 
0.0). 



gives rise to the parametrized counterpart of this model, 
namely the discrete Remoisseiret-Peyrard model. In this 
goal, we assume the usual parametrization of spin vari- 
ables using Euler angles i.e.: 



s„ — sin( 



I cos I 



= sm On smcpn, 



— COS ( 



where 0„ and (f>n are more precisely the local variables 
describing rotations of the n*'* spin. When these two 
angular variables assume arbitrary values, spins undergo 
precessional out-of-plane motions i.e. rotations out of the 
plane perpendicular to the chain axis. Nevertheless, we 
can confine spin motions in a fixed plane perpendicular to 
the chain axis by setting On = tt/2. With this constraint, 
the total Hamiltonian |^ gives rise to the following non- 
linear difference equations. 



sin{(j)n+i - 4>n) - sm(. 
1 sin <j)n 



[1 



2rcos ( 



= 0, 



2J 
/i(l -r2)2 



(16) 



(17) 



These equations describe the texture of the space 
spanned by the n^^ spin rotation in the presence of two 
nearest- neighbour spins. The energy density associate to 
this spatial texture is readily obtained as: 



COS[(f 
h{T) 



-1 - <t>n) + COs( 
1 - COS 4>n 



2 J 1 



2rcos ( 



(18) 



The continuum version of Hlfci|l . i.e. the Remoissenet- 
Peyrard equation, admits exact topological kink 
solutions. We can therefore anticipate predicting exact 
soliton solutions for the discrete counterpart. In this 
viewpoint, (|18|l is nothing else but the energy density 
needed for the kink soliton to keep its topological shape. 
We have simulated the discrete equation Hlt)|) 
as well as the energy density (|18ll . On fig- 
ures nil ca ini HI El cni El and m we piot 



results for r = 0.0, 0.2, 0.5, and 10, respectively. The 
most stricking feature in all the figures representing 
solutions is the well-formed kink shape. Moreover, the 
kink-lattice structure is well developped and apparently 
exists for a rather broad range of values of the model 
parameter. It is interesting to note that no kink solution 
exists for r = 1, which is predictable as there the 
"nonlinear force field" vanishes. To our knowledge, this 
is the only singular point of the model. In this respect, 
the usual restriction confining r between -1 and 1 is 
in deep contradiction with our own numerical findings. 
According to our simulations, the most remarkable 
fact at large values of r is the negative energy density. 
However, the kink-lattice structure still shows apparent 
stability as one can note on figure El Morever, both the 
dispersion relation and the anharmonic spin-wave 
modes of figure El showed that the ferromagnetic 
groundstate was not so drastically affected by large 
values of r. Instead, what we noted was some cycle 
of variation of physical parameters with increasing r, 
for instance a close look on figure reveals that the 
dispersion curve for r ~ 5 has a gap whereas the one 
corresponding to r ^ 1 or close is gapless. Similar 
behaviour also shows up in the variation of the shape of 
the ISL modes. Namely, figure [TUl suggests a waveshape 
not so different from those found for values of r less than 
1. The negative value of the energy density as given in 
figure El also relate to this particular variations of the 
model characteristic parameters at large r. However, 
for this last case, one will easily check that the negative 
energy density is compensated by increasing the ratio 
ho/ J- We can therefore expect similar corrections on the 
previously simulated quantities by judicious adjustments 
of this ratio. 



IV. CONCLUDING REMARKS. 

We have explored the static properties of a 
parametrized model as a probe of its ability to help 
describing the smooth variations of characteristic mag- 
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FIG. 12: Kink energy density for r = 0. 




FIG. 13: The kink lattice for r = 0.2. 

netic properties as it can occur in some ferromagnetic 
alloy series. In particular, we have in mind ferromag- 
netic systems with planar anisotropic structures but 
showing slight deviations from the conventional planar 
model, in which an applied magnetic field couples ex- 
actly to the spin component along the anisotropy axis. 
The parametrized model introduced as illustration dis- 
plays rather rich features to both viewpoints of its spin- 
wave phase and domain-wall phase. Recall that this 
parametrized model is actually an improved version of a 
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FIG. 14: Kink energy density r = 0.2. 
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FIG. 15: The kink lattice for r = 0.5. 
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FIG. 16: Kink energy density for r = 0.5. 



previously suggested spin model 0, 0| . While the cur- 
rent version is more rigourous i.e. more appropriate for 
ferromagnetic systems often involving single crystalline 
anisotropy, the original version which is more general can 
become quite powerful for systems in which both an easy- 
axis and an easy-plane anisotropics co-exist. Still, as far 
as we know, simultaneous crystalline anisotropics of this 
kind are more probable in ferrimagnets or antiferromag- 
nets but less in a true ferromagnetic spin system. 




FIG. 17: The kink lattice for r = 10. 
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FIG. 18: Kink energy density for r = 10. 
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